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SYSTEMS AND METHODS FOR CHARACTERIZING THE COVERAGE OF AD 

HOC SENSOR NETWORKS 

Statement of Government Funding 

[0001] The work described in this application was supported by funding 
from the Defense Advanced Research Projects Agency under DARPA contract 
number F30602-00-C-0139. The government may have certain rights to this 
invention. 

BACKGROUND OF THE INVENTION 

1. Field of Invention 

[0002] The systems and methods of this invention are directed to ad hoc 
sensor networks. 

2. Description of Related Art 

[0003] Ad hoc sensor networks have potential impact on many applications 
such as intelligent infrastructures, environmental monitoring, and battlefield 
monitoring. The characterization of a network's ability to perform a sensing task is a 
fundamental problem for the deployment and use of such systems. Traditional sensors 
are characterized by specifications such as range, resolution, and accuracy. These 
specifications are used to decide if a sensor is suitable for the task and system at hand. 
The simplicity and elegance of these specifications make the use of individual sensors 
a very easy task. Similar specifications for sensor networks currently do not exist and 
extensive testing is necessary to assess the capability of sensor networks and then their 
suitability for a particular task. 

[0004] One approach which tackles coverage problems of sensor networks is 
described by Meguerdichian et aL, "Coverage Problems in Wireless Ad Hoc Sensor 
Networks," IEEE Infocom 2001, vol. 3, pp. 1380-1387, April 2001. This publication 
uses the distance from the sensor nodes as the only metric that characterizes the 
sensing ability. It does not take into consideration any sensor characteristics. 
Accordingly, this results in the same coverage characterization independent of the 
type and capabilities of the individual sensors used in the network. Moreover, this 
approach computes a maximal breach path using a geometrical characterization that is 
not necessarily related to the sensing capabilities of the network. 
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[0005] U.S. Patent No. 6,212,471 to Stiles et al. discloses dynamic optical 
sensor employment for searching an area. U.S. Patent 6,297,763 to Thomson et al. 
discloses a method for optimizing the coverage area of a sensor. In both of these 
patents, coverage characterization of sensing systems is of great significance. 

SUMMARY OF THE INVENTION 

[0006] The systems and methods according to various exemplary 
embodiments of this invention use the concept of a sensing field for specifying the 
ability of a sensor network to measure a physical phenomenon of interest. In various 
exemplary embodiments of the systems and methods according to this invention, a 
sensing field is defined over the space in which the phenomenon of interest needs to 
be estimated. For example, if one is interested in measuring the temperature of an 
object, then the sensing field may be defined over a set of possible temperature values. 
If, for example, one is interested in some phenomenon in a physical location domain, 
the sensing field may be defined with respect to that domain. The value of the sensing 
field at any given point reflects the ability of the sensor network to estimate the 
physical phenomenon of interest at that point. 

[0007] Various exemplary embodiments of the systems and methods 
according to this invention use a statistical method that determines a sensing field 
based on sensor layouts and the sensing models, thereby defining well-monitored 
regions and sensor holes in the sensor network. This information is used to 
characterize the quality of service that the sensor network provides for different 
applications. A graphical user interface may be used to display this information to the 
user for monitoring and health management of the sensor network. 

[0008] One object of the systems and methods according to this invention is 
to specify coverage of a sensor network using a sensing field. 

[0009] Another object of the systems and methods according to this 
invention is for the sensing field to characterize the ability of the sensor network to 
measure a physical phenomenon. 

[0010] Various exemplary embodiments of the systems and methods 
according to this invention display the sensing field, for example, using a graphical 
user interface, thereby providing a simple and intuitive way to monitor the coverage of 
the sensor network. 
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[0011] Various exemplary embodiments of the systems and methods 
according to this invention include specification of the coverage using a sensing field 
that takes into consideration the sensor layout and the sensing models. Given the 
sensor network layout and sensing models, the resulting sensing field coverage 
obtained by the exemplary embodiments of the systems and methods of the invention, 
accurately represents the ability of the sensor network to measure physical 
phenomena. 

[0012] The sensing field aspect of the various exemplary embodiments of 
the systems and methods according to this invention may be used to solve various 
problems in sensor networks, such as, for example, sensor deployment, sensor 
localization, and sensor health management. The sensing field may also be used from 
an adversarial point of view, to find weak points and escape detection. 

[0013] Other features and advantages of this invention will become apparent 
from the following detailed description of exemplary embodiments. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0014] Various exemplary embodiments of this invention will be described 
in detailed with reference to the following drawings, wherein: 

[0015] Fig 1 is an example of a sensor observational model that plots 
target location distribution based on a single amplitude measurement: 

[0016] Fig. 2 is an example of a sensor observational model that plots the 
target location distribution based on a direction-of-arrival measurement. 

[0017] Fig. 3 is a plot of overall target location distribution according to an 
exemplary embodiment according to the systems and methods of this invention; 

[0018] Fig. 4 is a plot of a best Gaussian fit to the overall target location as 
shown in Fig. 3; 

[0019] Figs. 5 is an example of a plot of the sensing uncertainty, measured 
as the determinant of covariance matrix, as a function of target distance according to 
an exemplary embodiment of the systems and methods of this invention; 

[0020] Fig. 6 is an example of a plot of the goodness of fit, measured as the 
Kullback-Leibler distance as a function of target distance, according to an exemplary 
embodiment of the systems and methods of this invention; 
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[0021] Fig. 7 is a flowchart outlining one exemplary embodiment of a 
method for choosing a network arrangement according to the systems and methods of 
this invention; 

[0022] Fig. 8 is a block diagram of an exemplary embodiment of a sensor 
network system with combined sensing ability characterization, visualization, and 
application selection according to this invention. 

[0023] Fig. 9 is an example of a sensor network which has heterogeneous 
sensor types; 

[0024] Fig. 10 is an exemplary plot showing the localization uncertainty 
using all sensors in the sensor network of Fig. 9; and 

[0025] Fig. 1 1 is an exemplary plot showing the localization uncertainty 
using only four nearest sensors in the sensor network of Fig. 9. 

DETAILED DESCRIPTION OF EXEMPLARY EMBODIMENTS 

[0026] In one exemplary embodiment of event localization problem solving, 
according to the systems and methods of this invention, a network with heterogeneous 
sensors, range, and angle is considered, and a sensing field is determined that specifies 
the ability of the sensor network to estimate the position at every location. Various 
exemplary embodiments of the systems and methods according to this invention make 
such a determination based on the uncertainty of the estimation and the quality of the 
sensing model. The notation of estimation uncertainty and sensing quality is 
universal over various sensor network applications; therefore, the sensing field and the 
determination described herein are not limited to any particular application, but rather 
apply generally to a wide variety of sensor network applications. 

[0027] In exemplary embodiments of event localization, such as, for 
example, target localization, according to the systems and methods of the invention, it 
is reasonable to assume that the target is present in a location domain. The domain 
may be, for example, a one dimensional road, a two-dimensional terrain, or a three- 
dimensional space. One goal of localization is to locate the target in the domain with 
some accuracy based on the observed data. Sensors are used to gather information 
about the target from various locations. By exploiting the spatial diversity of the 
sensors, one can achieve good localization performance, which would otherwise be 
unachievable using any single sensor. 
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[0028] Various exemplary embodiments of the systems and methods 
according to this invention determine in which regions an event, such as, for example, 
a target, can be easily detected, and in which regions the event, e.g., target, can 
possibly hide, even though the target has never been present in the location domain, 
based on the sensor layouts and the sensing models. Information regarding the sensor 
layouts and the sensing models is very useful for randomly deployed sensor networks. 
Once the sensor locations are calibrated, a user may obtain a good understanding of 
the strengths and weaknesses of the sensor network's sensing capabilities, and thus 
may adapt the sensor network accordingly. Various exemplary embodiments of the 
systems and methods according to this invention also permit the user to adjust the 
sensor layout to provide even coverage across a domain or entire field, or to focus 
coverage to a specific region of interest, depending upon practical situations. 

[0029] Sensor networks usually contain a large number of sensors spread 
over a region. For each sensor, a sensing field may be defined. To define a statistical 
tracking system that characterizes the pointwise ability of a network of sensors to 
measure a physical phenomenon of interest, the following two factors may be 
considered. The first factor is that the larger the uncertainty of the physical 
phenomenon, such as the location of the target in localization applications, the weaker 
the sensing ability. The second factor is the degree to which the actual sensing model 
fits the statistical assumptions. Various exemplary embodiments of the systems and 
methods according to this invention determine and illustrate these two factors with 
examples. The definition and determination of the sensing ability field according to 
exemplary embodiments of the systems and methods of this invention are applicable 
to, but not restricted to, the examples discussed below. 

[0030] To determine the first factor, the systems and methods according to 
this invention in one exemplary embodiment use Lo to denote the true target location, 
and L to denote the sensor measurement of target location. Uncertainty reflects the 
compactness of the conditional density p(L / Lo). For a one-dimensional case, the 
uncertainty can be measured by the variance of p(L / Lo). For higher dimensional 
cases, the determinant or the trace of the covariance matrix may be used. The former 
measures the volume of the probability mass, whereas the latter measures the 
circumference of the probability mass. Depending upon the application, a user may 
choose the appropriate measure. In general, the higher the uncertainty, the more 
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difficult it is to estimate L 0 from L. [0031] The second factor may be expressed by 
the goodness of fit of statistical assumptions. The statistical system according to one 
exemplary embodiment makes exemplary simplifying assumptions about physical 
reality. For example, with respect to analysis and computation, sensing models often 
are provided in parametric forms. For tracking, a Gaussian assumption is common. 
For example, Kalman filtering assumes Gaussian observation modes. Linear 
minimum mean-squared error (LMMSE) estimators are optimal for Gaussian cases, 
and serve well as an approximation for non-Gaussian distributions. In reality, the 
actual sensing model may deviate from the assumptions. To measure the difference 
between an assumed model q(L /l 0 ) and the actual model p(L / Lo), the information- 
theoretic Kullback-Leibler difference may be used. The Kullback-Liebler difference, 
is defined in equation 1, below. 

[0032] KL(p, q) = j log(p/q) dp (1) 
[0033] The Kullback-Leibler divergence is known to those skilled in the 
mathematical arts and is described, for example, by S. Kullback, Information Theory 
and Statistics, John Wiley & Sons, New York, 1959. The Kullback-Liebler 
divergence measures the number of bits to describe the discrepancy between p and q. 
Alternatively, the Hellinger divergence, which is specified by equation 2, below, and 
which is similar to a mean-squared error (MSE) measure, may be used. 

[0034] H(p, q)= J[sqrt(p)- sqrt(<?)] 2 dp (2) 
[0035] These metrics tell how well the assumption q (L / Lo) fits the actual 
sensing physics. A good fit would justify the use of a tracking algorithm built on the 
assumption q. A bad fit would suggest otherwise. 

[0036] Based on a consideration of these two factors, a simplified exemplary 
embodiment is described of the systems and methods of this invention is described 
below, which assumes that a sensor is located at the origin (0,0) without loss of 
generality, and that the sensor is microphone based, with a measurement of the 
amplitude Z and a measurement of the direction-of-arrival (DOA) measurement data. 
[0037] The physical observation model consists of the following: 
[0038] First, conditioned on a true location Lo, the amplitude measurement Z 
and the estimated direction of arrival 0 are independent. 

[0039] Second, given Lo, Z= A/r(Lo) + «, where r(L 0 ) is the distance from Lo 
to the sensor location (the origin), A is the sound source amplitude, and n is the 
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measurement noise, assuming that the noise is independently and identically 
distributed (iid) with a Gaussian distribution with zero mean, and variance a n 2 . From 
this amplitude measurement, the sensor infers the target distance as r = A/Z. 

[0040] Third, given L 0 , 9 = 0o+w, where 9 0 = atan2(L 0 , y , L 0fX ) is the angle 
from the sensor to the true target location, and w is measurement noise, assuming that 
the noise is independently and identically distributed (iid) with a Gaussian 
distribution with zero mean, and variance a w 2 . 

[0041] The target location may be determined from the range and direction 
measurements as L x = r cos(0) and L y = r sin(O). From the sensing field model, the 
conditional probability may be derived as: 

p(L Il 0 ) =p(L x ,L y lU) =p(r,6 /L 0 )r = V {r \ L 0 ) • ^L x 2 +L y 2 

[0042] Figs. 1 and 2 plot e one example of a distribution of p(r / L 0 ) and p(0 
I Lo) as a function of L x and L y , for Lo =(30, 0), A = 80, a n = 0.5, and <r w = 5 degrees. 
The overall distribution p(L / L 0 ) is plotted in Fig. 3. Note that it appears to be a 
comet around the true location Lo. Fig. 4 is the Gaussian fit of the aforementioned 
overall distribution. Compared to the actual distribution, the Gaussian fit is 
symmetric and centered at L 0 . In this case, the Gaussian fit approximates the original 
distribution p(L I L 0 ) reasonably well. 

[0043] Given this example of the sensor model, the sensing field may be 
determined. The sensing field should be symmetric around the sensor location. 
Accordingly, only the target locations L 0 = (lo, 0) for L 0 in the positive half axis need 
be determined. For illustration purposes, the uncertainty may be measured as the 
determinant of the covariant matrix. Fig. 5 illustrates that the uncertainty increases 
with the radius Lo. For this sensing model, when the target is far away, the signal is 
buried in noise. Therefore, the location measurement is less accurate. 

[0044] The goodness of fit may be measured as the Kullback-Leibler 
divergence, as illustrated in Fig. 6. This divergence is small for a mid-range of L 0 , 
and large outside the mid-range. For a small L 0 , since the target is very close by, the 
signal-to-noise ratio (SNR) is high. Therefore, p(r / L 0 ) is essentially a very thin ring, 
and p(L / L 0 ) is a thin arc which cannot be well approximated by a Gaussian fit. For 
very large L 0 , the tail of the comet becomes longer, and the Gaussian fit is not 
accurate. 
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[0045] The sensing ability characterization for single sensors may be 
combined to generate a sensor ability field for a sensor network which includes 
multiple sensors. The Kullback-Leibler divergence field of each sensor may be 
superimposed to measure the goodness of fit. The uncertainty of target location with 
respect to combination of sensors is measured as the determinant of the covariance 
matrix, and the sensing quality is inversely proportional to the uncertainty. Hence, the 
target location uncertainty is measured as l/det(cov) These approaches may be 
justified from a Bayesian estimation point of view. The regions with low uncertainty 
and low divergence values are considered to be well monitored by the sensor network. 
The regions with high uncertainty and high divergence values are considered to be 
sensor holes. In this manner, the ability of the sensor network may be accurately 
characterized. 

[0046] The systems and methods according to this invention are not 
restricted to target localization issues. Various exemplary embodiments according to 
the systems and methods of this invention which characterize a sensor network's 
ability by measuring the uncertainty of the estimate and the goodness of fit of 
statistical assumptions, although fundamental in sensing applications, may also be 
applied to very different scenarios such as temperature estimation or chemical 
monitoring. Further, various exemplary embodiments of the systems and methods 
according to this invention may be applied to not only fixed sensors, but to mobile 
sensors as well. 

[0047] For mobile sensors, sensor locations may be acquired from some 
location service, for example, from the readings of GPS devices. The location and 
sensing model information can be periodically aggregated to some central 
computation center, and by the same method, the sensing field can be computed. 
Assuming the sensor nodes move in a continuous fashion, one may design incremental 
updates to save computation in the updating of the sensing field. 

[0048] Fig. 7 is a flowchart which outlines one exemplary embodiment of a 
method according to this invention for choosing a network arrangement based on a 
sensing field quality characterization. Control starts in step SI 000 and proceeds to 
step SI 010, where a sensing field is defined for each sensor node of the field. Control 
then proceeds to step SI 020, where a determination is made of well monitored regions 
and sensor holes in the sensing field. 
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[0049] Control then proceeds to step SI 030, where the quality of the 
network for various applications is characterized. 

[0050] Control then proceeds to step SI 040, where the quality 
characterization from the previous step is visually displayed, for example, on a 
graphical user interface. Control then proceeds to step SI 050, where a network 
arrangement based on the sensing field quality characterization previously determined 
is chosen. Control then proceeds to step SI 060, where the process ends. 

[0051] Various exemplary embodiments of this invention may be 
implemented on a programmed general purpose computer. However, various 
exemplary embodiments can also be implemented on a special purpose computer, a 
programmed microprocessor or microcontroller and peripheral integrated circuit 
elements, an ASIC or other integrated circuit, a digital signal processor, a hardwired 
electronic or logic circuit such as a discrete element circuit, a programmable logic 
device such as a PLD, PLA, FPGA or PAL, or the like. In general, any device, 
capable of implementing a finite state machine that is in turn capable of implementing 
the flowcharts shown in FIG. 7, can be used to implement systems and methods of this 
invention. 

[0052] Moreover, various exemplary embodiments can be implemented as 
software executing on a programmed general purpose computer, a special purpose 
computer, a microprocessor or the like. In this case, the systems and methods of this 
invention can be implemented as a routine embedded in a printer control system or 
controller or the like. That is, the systems and methods of this invention can be 
implemented by physically incorporating them into a software and/or hardware 
system, such as the hardware and software systems of a personal computer. 

[0053] Fig. 8 shows an exemplary embodiment of a system according to 
this invention. A sensor network 100 is placed in communication with a link or bus 
200 via a link 101 . Further, a sensing application, circuit or routine 300, a 
determination application, circuit or routine 400, a characterization application, circuit 
or routine 500 and a selection application, circuit or routine 600 are placed in 
communication with the bus 200 via links 301, 401, 501 and 601, respectively. 

[0054] As discussed above, each sensor node is relatively self sustaining in 
the sense that each sensor node has its own power supply. Each sensor node also has 
a central processing unit, memory and communication (transmission and reception) 
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capabilities. Each sensor node may communicate with the sensing application, circuit 
or routine 300, the determination application, circuit or routine 400, the 
characterization application, circuit or routine 500 and the selection application, 
circuit or routine 600 via the bus 200. These applications, circuits or routines are 
capable of carrying out the functions described above with respect to Fig. 7. 

[0055] A display device 700 and an input/output device 800 may also be 
placed in communication with the bus 200 via links 701 and 801, respectively. The 
display device 700 may be used to display the sensor field quality. The input/output 
device 800 allows a user to interact with the various applications, circuits and 
routines. Further, the input/output device 800 may allow a user to interact, either 
directly or indirectly, with the sensor network 100. 

[0056] It should be appreciated that each of the applications, circuits or 
routines shown in Fig. 8 can be implemented as portions of a suitably programmed 
general purpose computer. Alternatively, each of the applications, circuits or routines 
shown in Fig. 8 can be implemented as physically distinct hardware circuits within an 
ASIC, or using a FPGA, a PDL, a PLA, or a PAL, or using discrete logic elements or 
discrete circuit elements. The particular form each of the applications, circuits or 
routines shown in Fig. 8 will take is a design choice and will be predictable to those 
skilled in the art. 

[0057] Further, it should be appreciated that each of the links shown in 
Fig. 8 can be implemented as any known or hereafter developed technology or 
technique that allows communication between components. In particular, the links 
may be wired or wireless. 

[0058] Fig. 9 shows a two-dimensional network sensor layout. The points 
marked with "+" denote amplitude sensors, such as, for example, acoustic sensors, and 
assume that the amplitude A equals 40 and the Gaussian amplitude noise of standard 
deviation a k , a m P is 0.5. The points marked with squares denote direction-of-arrival 
(DOA) sensors, and assume that the noise-to-angle measurements have Ok,DOA is three 
degrees. 

[0059] Fig. 10 shows an uncertainty field for the exemplary sensor network 
of Fig. 9. The uncertainty is measured as the determinant of the Fisher information 
matrix using all sensors in the sensor field. Fig. 10 is plotted in log scale. 
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[0060] Fig. 1 1 shows another uncertainty field using a determination of an 
instantiation of sensing field based on the Cramer-Rao bound (CRB) of estimation 
accuracy, or, equivalently, the Fisher information, in sensor networks. Fig. 1 1 
involves the determinant of the distributed Fisher information. Details of the 
determination may be found in Liu, Juan, et al., "Sensing Field: Coverage 
Characterization in Distributed Sensor Networks", Proceedings of ICASSP, May 
2003, the entirety of which is hereby incorporated by reference. Fig. 1 1 is plotted in 
log scale. 

[0061] Sensing field determinations for characterizing the pointwise ability 
of a sensor network to measure and/or otherwise determine a physical phenomenon 
of interest may be used as a tool to evaluate the sensor network, and provides a 
reference to sensor deployment as well as estimation and routing algorithms in sensor 
networks. 

[0062] For example, for a randomly deployed sensor network, once the 
sensor locations/nodes are calibrated, a user can obtain a good sense of the strength 
and weakness of the sensor network's sensing capabilities, and algorithms can be 
adapted accordingly. Furthermore, the system's awareness of a target fast approaching 
a hole in the sensor network may trigger algorithms to trigger sensors further away, 
call for more extensive sensor collaboration, and planning for limited resources with 
sufficient look-ahead capabilities. 

[0063] As another example, in a mobile sensor network, where adjusting 
sensor layout is important, individual sensors can be adjusted to provide relatively 
even sensor coverage across the field, or to focus coverage on specific regions of 
interest to satisfy user requirements. 

[0064] While this invention has been described in conjunction with the 
exemplary embodiments outlined above, it is evident that many alternatives, 
modifications and variations will be apparent to those skilled in the art. Accordingly, 
the exemplary embodiments of this invention, as set forth above, are intended to be 
illustrative and not limiting. Various changes may be made without departing from 
the spirit and scope of this invention. 



